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Abstract. We discuss two modifications of domain-wall fermions, aimed to 
reduce the chiral-symmetry violations presently encountered in numerical 
simulations 1 . 



1. INTRODUCTION 

Domain-wall fermions (DWF) live on a five-dimensional lattice [1, 2, 3]. The 
fifth coordinate, usually denoted s, is introduced because chiral-symmetry 
violations (xSV) decrease with increasing size of the fifth dimension, N s . 

DWF are rapidly becoming a standard method for lattice QCD sim- 
ulations (for a review see ref. [4]). While showing agreement with other 
methods for quantities such as the strange-quark mass and Bk [5], a re- 
cently reported first result for e'/e [6] signals that DWF may be employed 
to compute the Standard-Model predictions in cases where the traditional 
(Wilson, staggered) fermion methods have so far been unsuccessful. DWF 
have also been employed for the study of thermodynamics, topology and 
the chiral condensate [7]. We note that the four-dimensional overlap-Dirac 
operator [2, 8] arises in a special limit of DWF [9], hence results pertaining 
to this operator have relevance for DWF too. 

When the inverse lattice spacing is a~ l ~ 2 GeV (quenched (3 ~ 6.0) 
%SV are found to be at the few percent level for N s ~ 15 — 20 for most 
quantities [10, 5, 6, 11]. However, at stronger coupling, a" 1 ~ 1 GeV, the 
violations are much bigger [7, 11], and sometimes even N s ~ 50 does not 
seem to be enough. Since the cost of simulations grows linearly with N s , 
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an obvious question is whether fiddling with the DWF action may allow us 
to decrease N s while maintaining (or even improving) the quality of chiral 
symmetry. 

In this report we describe two modifications of DWF, aimed to reduce 
the presently encountered xSV. The first method, discussed in Sect. 2, is 
applicable in dynamical simulations [12]. This method is based on non- 
perturbative considerations. The second method (Sect. 3) is based on per- 
turbative considerations, and is applicable in both quenched and dynamical 
simulations [13]. 



2. A MODIFIED PSEUDO-FERMION ACTION 

The DWF determinant can be written as 

det (D F ) = (hf) Ns x (finite factor) . (1) 

Here Dp is the five-dimensional DWF matrix (we use the notation of 
ref. [12] with lattice spacings a = 05 = 1). The "finite factor" has a conver- 
gent N s — > 00 limit. It contains the physics of the quark field, whose two 
chiral components live on opposite four-dimensional boundaries of the lat- 
tice. The bulk factor (/j-f) Ns arises because of the presence of extra N s — 1 
four-dimensional fermions with masses of the order of the cutoff. This unde- 
sirable factor is canceled by introducing opposite-statistics Pseudo Fermion 
(PF) fields [2, 14]. In the so-called second order formulation, the PF action 
is 

S pi = Qxs ( D pf D if)xs,ys> <t>ys> , (2) 

xs,ys' 

where the PF lattice has only jV p f = N s /2 sites in the fifth direction. 
The standard choice is D p f = Dp(m = 1). The parameter m controls 
the coupling between the two boundaries of the five-dimensional lattice 
(hence also between the left and right components of the quark field), and 
is proportional to the bare quark mass for m <C 1. The choice m = 1 
corresponds to anti-periodic boundary conditions, and implies the absence 
of any light PF states. 

Aiming to reduce %SV we consider a modified PF action 

xs,ys' 

+ cV x (3) 

xs,ys' 

where V x = E M<J y Retr (^ ~ U x ^U x+ ^ u Ul + ~^Ul u ) is the plaquette-action 
density at the point x, and c is a coupling constant. 

In Subsect. 2.1 below we discuss what we believe to be the primary 
source of xSV at strong coupling. In Subsect. 2.2 we explain how the mod- 
ified PF action should help reducing xSV. Our considerations also suggest 
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how to look for the parameters' range where the new PF action may be ef- 
fective. The new term in eq. (3) does not change the sparseness of the (2nd 
order) PF matrix. Therefore the extra cost of simulating it is hopefully not 
large. (One should give thought to finding a good pre-conditioning.) 

2.1. THE PROBLEM 

It was first observed by Narayanan and Neuberger [2] that the DWF de- 
terminant can be compactly expressed in terms of a (second-quantized) 
transfer matrix T that describes single-site hopping in the unphysical s- 
direction. The T-formulation is especially powerful for DWF, because the 
same gauge field couples to the fermions on every s-slice, and so T is inde- 
pendent of s. 

As expected from their good chiral properties, DWF have (on a finite 
lattice, almost) exact zero modes in an instanton background. As a topo- 
logically non-trivial gauge field is gradually turned on, at some point level- 
crossing occurs in the spectrum of the background-field dependent hamil- 
tonian H = logT, and an initially-excited state becomes the new ground 
state [2, 15, 16]. Beyond the crossing point the axial U(l) charge is not 
conserved, as required by the anomaly. At the crossing point, one of the 
eigenvalues of T, the corresponding first-quantized transfer matrix, is equal 
to one. Two comments are in order. 1) An eigenvector of T with eigenvalue 
one exists if and only if there exists a zero eigenvector of the hermitian 
Wilson-Dirac operator (this extends to a mapping between all eigenvectors 
of T and of the hermitian Wilson-Dirac operator [12]). 2) Since the mass 
of the Wilson-Dirac operator is supercritical, and is far from the critical 
mass M c , the above zero eigenvectors have nothing to do with continuum 
physics. They are highly localized states whose support typically extends 
over only one or two lattice spacings [17]. Moreover, one can rigorously 
prove the absence of such zero eigenvectors if the plaquette-action density 
is everywhere smaller than a certain constant [18]. 

A near-unity eigenvalue (NUEV) of T implies little or no suppression 
of s-correlations for quark states that have a finite overlap with the cor- 
responding (four-dimensional) eigenvector. Strong correlations across the 
s-direction also mean strong lattice-artifact %SV for (non-singlet) axial 
currents [3]. As follows from the picture of topology changing developed 
in ref. [2], NUEVs cannot be completely avoided. This said, the spectral 
density of NUEVs is a very sensitive function of the lattice parameters. 
The dependence of the spectral density on [3 was investigated in ref. [19] 
for quenched ensembles. The results (see in particular fig. 5 therein) confirm 
that NUEVs are extremely rare for j3 > 6.0. Practically, this implies an ex- 
ponential fall-off of all correlations in the s-direction. However, at stronger 
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coupling there is a dramatic rise in the density of NUEVs. In this case xSV 
decay very slowly and are significant even for very large N s [11]. In fact, it 
was recently proved within the strong-coupling expansion that explicit xSV 
persists in the limit N s — > oo, and that the symmetry group of the effective 
hamiltonian is identical to that of naive fermions with a mass term [20] . 

What is the role of NUEVs in a dynamical simulation? As we explain 
in Subsect. 2.2 below, to every NUEV of T corresponds a family of small 
eigenvalues of the DWF operator. If "small" means, say, O(Aqcd), the 
number of small eigenvalues of Dp is 0(N S Aqcd)- Why, then, are the 
relevant "bad" gauge- field configurations not suppressed for N s S> 1? The 
same set of small eigenvalues occurs also in the PF determinant which 
sits in the denominator, and the two sets simply cancel each other. This 
is a characteristic feature of small eigenvalues of Dp which are related to 
NUEVs of T. In contrast, the number of DWF (near) zero modes found in 
an instanton background is as required by the anomaly (not proportional 
to iV s ) and with no matching PF zero modes. 

One can say that "bad" gauge-field configurations are not suppressed 
because the PF determinant follows the DWF one too closely! However, 
there is no need to maintain this tight relation far from the continuum 
limit. The purpose of the new PF action eq. (3) is to break this relation in 
a controlled way, such that configurations with NUEVs will be suppressed. 



2.2. A POSSIBLE SOLUTION 

We now turn to the analysis of the modified PF action. (The PF action (3) 
was first introduced in ref. [12]. The below discussion is largely new. While 
we try to make this report self-contained, some familiarity with ref. [12] 
should help in keeping track of what follows.) The crucial step is to rewrite 
the PF determinant as a product, det Npf (B) det(f2(c)), where 

&{c)xs,ys' = -U s ,s' &x,y + S s>s ' ^a{c) x ^ + boundary terms , (4) 

Ot(c) = 75£T 1/2 [Dpf^Jf + cV] ^B~ l/2 . (5) 

Here D is the massless Wilson-Dirac operator, V is the diagonal matrix with 
entries 5 x ,yP x , and B = 1 — M + W where W is the Wilson term and M is 
the five-dimensional mass, or "domain-wall height". (Both M and W are 
assumed to be positive; because of the relative minus sign this corresponds 
to supercritical Wilson fermions.) The free nearest-neighbor laplacian for 
the s-coordinate is D s , s ' = $s,s'+i + $s,s'-i ~ 2$s,s'- The boundary terms in 
eq. (4) are zero unless s or s' are equal to 1 or AT pf . We will be interested in 
computing how the Boltzmann weight varies as a function of c, i.e. in the 
ratio 

exp(-ft(c)) = det(n(0))/det(fi(c)) . (6) 
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Since the (PF) action is local, the free energy is dominated by a bulk term 
that grows linearly with N p {. We will only be interested in the way the 
PF bulk factor changes. Writing 11(c) = N p{ K^c) + K (c) + 0(1/N p{ ) we 
observe that that 1Z\(c) is independent of the boundary terms in eq. (4) 
(this can be proved using the determinant formulae of ref. [12]). Denote 
by O'(c) the operator that coincides with 0(c), except that the boundary 
terms are absent (i.e. O'(c) lives on a periodic lattice in the s-direction). 
The eigenvectors of Jl'(c) are given by e %SVh ^f\(x) where e lsp -> is a one- 
dimensional plane wave and ^\(x) is an eigenvector of 04(c) with eigenvalue 
A(c) (04(c) is positive, hence A(c) > always). 

The eigenvalues of O'(c) are given by A(c) + 2(1 — cos(ps)). Consider 
first the situation at c = 0. If the transfer matrix has a unity eigenvalue, 
^4(0) has a zero eigenvector ipo, and therefore O'(0) has a gap-less family of 
eigenvalues given by 2(1 — cos(ps)). This explains the statement made in the 
previous subsection about the number of small DWF and PF eigenvalues 
for "bad" configurations. Projected [12] onto the four-dimensional eigenvec- 
tor ipo, the DWF propagator G(s,s') = J2 P5 e ips ^ s ~ s ) (1 — cos^)) -1 does 
not fall exponentially with the separation \s — s'\, leading to unsuppressed 
correlations across the s-direction. 

In computing 1Z\ (c) we will simplify things by replacing the lattice cutoff 
for the s-direction with a sharp momentum cutoff. Since the modification 
of the PF action is intended to be small, it should affect mainly the small 
eigenvalues, and changing the cutoff scheme should be harmless. As a second 
simplification we will use first-order perturbation theory to determine how 
the eigenvalues of 04(c) vary. One has A(c) ps Aq + cAi where 



(note that Ai > 0). Replacing the fifth-coordinate momentum sum by an 
integral (which is allowed for large N p f) yields the following result 



We have sent the momentum cutoff to infinity. The result is finite, and this 
justifies our treatment of the s-coordinate's ultra-violet cutoff. 

We now wish to distinguish two, qualitatively different, cases depending 
on the gauge-field background. First assume that there is a gap in the 
spectrum of O 4 (0), namely Ao > A m i n > for all Ao where by assumption 
Amin = 0(1). Taking c <C A m i n (note that Ai = 0(1) by eq. (7)) and 





U ^ c ) = E/ S ( 1o §(P5 + Aq + cAi) - log(^ + Ao)) 




(8) 
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expanding the first term in eq. (8) we have 

^( c ) = c E^r- ( 9 ) 

What is the physical significance of this result? We can compute lZ\{c) 
also by varying the PF partition function (defined with periodic boundary 
conditions) with respect to c. This yields 

n 1 (c) = c J2r x (4s^s) , (io) 

X 

(for any fixed s). One can show that eqs. (9) and (10) agree for c <C A m i n . In- 
voking translation invariance i.e. replacing ((frxs^xs) by its four-dimensional 
average V^~ l J^xi&xs&xs) tells us that the effect of the modified PF action 
is an additive renormalization of the inverse bare coupling 

p +cNp{V -^(4s<t>xs) ■ (ii) 

X 

An estimate of the expectation value in eq. (11) may be obtained using 
standard weak-coupling perturbation theory. 

The second possibility is that the spectrum of ^4(0) is gap-less. For 
the low eigenvalues of ^4(0) we can no longer use eq. (9). Instead, for any 

Ao <C c one has 

VAo + cAi - v/A^- VcAi- (12) 

We see that, in the absence of a gap, lZ\{c) is not a linear function of c. 
This happens because of the singularity that det(f2(0)) develops when an 
eigenvalue Ao tends to zero. For small c, the r.h.s. of eq. (12) is much big- 
ger than c Ai/VAmin- Thus, as desired, no-gap configurations are relatively 
suppressed! 

In practice one can try to implement this new PF action along the fol- 
lowing lines. Suppose that we wish to reduce xSV in a dynamical simulation 
with given (3 = (3q. We start instead with a smaller value, say in the range 
(3 ~ 0.8/?o - 0.9/3o, and increase c until the measured lattice spacing agrees 
with the one previously obtained for j3 = (3q . At this point one should have 
a coupling constant c of order 0(1/N S ). (When varying N s it is advisable 
to keep the product cN s invariant by taking c = c\/N s with fixed c±.) For 
configurations with a gap the Boltzmann weight should now be roughly 
the same as before, because the reduction in the explicit (3 is matched by 
the effective increase in (3 coming from the PF action. On the other hand, 
no-gap configurations should now have a much smaller Boltzmann weight, 
since the suppression provided by the new PF action wins over the effect 
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of working at smaller [3. Thus, if there is a "window" where all the sim- 
plifications made above approximately hold, this method should lead to 
significantly reduced xSV. 

2.3. ADDENDUM 

We close this section with two comments. The favorable value of M is 
known to be M ~ 1.8. This implies that B = 1 — M + W is not necessar- 
ily a positive operator and, consequently, that the transfer matrix is not 
necessarily positive either. In practice, it turns out that the positivity of 
B is preserved. For example, for an ensemble of quenched (3 = 5.7 con- 
figurations, the spectrum of B was found to be bounded away from zero 
while M < 2, so the Wilson term has an effective lower bound of around 
one [21]. A similar shift is found for M c , the critical Wilson-fermion mass, 
which is also compatible with the mean-field and one- loop predictions [22] . 
We thus have a nicely consistent picture of several quantities all having the 
same origin - the rise in the effective lower bound on W away from the 
continuum limit. 

The second comment corrects a misleading statement made in ref. [12]. 
There, concern was expressed that a slightly different transfer matrix de- 
noted T could have complex eigenvalues. This actually cannot happen [23] 
because the transfer matrices T and T have the same spectrum. When B 
is negative, T may indeed have negative eigenvalues, but cannot have com- 
plex eigenvalues since its definition involves only B and its inverse, but not 
y/~B. (Having negative eigenvalues for T could also be annoying, since when 
log(T) does not qualify as a hamiltonian there could appear new unwanted 
effects; this unpleasant situation, too, seem to have little relevance because 
of the reasons explained in the previous paragraph.) 

3. THE FOUR-DIMENSIONAL PART OF THE ACTION 

In this section we propose a different modification of the DWF action, 
which may be used in both quenched and dynamical simulations (possibly 
in conjunction with the modified PF action in the latter case; the full details 
will be reported separately [13]). We begin by noting that, for a semi-infinite 
s-coordinate, the quark state with zero four-momentum attached to the 
s = 1 boundary falls like x( s ) ^ 1 1 — iV^| s . (This should not be confused with 
the "quark field" introduced in ref. [3]. The latter is an interpolating field 
for quark states which by definition is restricted to the boundary layers.) 
On a finite lattice, the mixing of the quark states from the two boundaries 
results in a residual bare quark mass m rcs = 0((1 — M) Ns ), even if the 
explicit bare mass m is zero [3, 24]. 
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In a simulation, the measured pion-mass squared will be proportional 
to the sum m + m res (as long as both m and m res are small) but m res 
is much larger than suggested by the tree- level result. Assuming that low 
momentum quark states behave in simulations like x( s ) ~ Qefi> one expects 
m res oc q^g . Thus, we can extract q e fi from the iV s -dependence of the pion- 
mass squared (or of the anomalous term in the relevant Ward identity). 
Results at quenched (3 = 6 [10] suggest a value q c s ~ 0.8. This value is very 
large, if we remember that q e s = 1 means no exponential fall-off. 

In the free DWF case, one has qo = |1 — M\ — ► for M — ► 1. Also in a 
mean- field approximation [22] one finds a similar result g m f = \\+8M — M\. 
If mean-field was a good description, letting M — ► 1 + <5M would give 
(7 e ff <C 1 in simulations. Since in practice q e g is much closer to one than it 
is to zero, we conclude that the effective value of SM varies considerably 
over different configurations and, in fact, over different spacetime regions 
of the same configuration. 

We see that tree-level or mean-field approximations fail to describe a 
key feature of the quark's wave function. One reason for this failure, namely 
NUEVs, was already discussed in the previous section. However, we believe 
that NUEVs cannot be the whole story. A large density of NUEVs implies 
no exponential fall-off at all, and a non-zero residual mass even for N s > 50. 
On the other hand, for quenched (3 > 6.0 the exponential fall-off seems to be 
there, albeit with a pretty large q e g. As we will now explain, we believe that 
the large value of q e g is essentially generated by perturbative fluctuations of 
the gauge field (this is also supported by the results of ref. [24, 16]). Below, 
we use the one-loop expression for the DWF self-energy derived in ref. [22] 
(for one-loop results for composite operators see ref. [25]). 

3.1. THE PERTURBATIVE WAVE FUNCTION 

Let us assume m = and N s 3> 1, so that the tree-level value of m res is 
negligible. For small four-momentum the free DWF propagator in the 
vicinity of the s = 1 boundary is 

Go(p„; s > s ') = Xo(s) - xo(s') Pl + Reg. (13) 

where Xo( s ) = |1 — M\ s is again the tree-level quark's wave function, Pl = 
5(1 — 75) and "Reg" stands for a regular function of p^. 

At the one-loop level, the leading quantum effect is the additive renor- 
malization M — > M + SM mentioned above [22]. This effect comes from a 
tadpole self-energy diagram, and must be treated non-perturbatively ("tad- 
pole improvement"). Having done so, we obtain the resummed one- loop 
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propagator [26] 

8 , s') = xi(s) ^ xi(s') + Reg. (14) 

where Z = 1 + 0(g 2 ) is a standard wave function renormalization factor. 
At the one-loop level we find the wave function 

Xi(s) = \l + SM-M\ s + g 2 S X i(s), (15) 

<5xi(s)octrP L S(^ = 0;s / = l,s) oc s~ 2 , (16) 

where S(/c;s',s) is the DWF self-energy coming from the "setting sun" 
diagram [22] (tr Pl s', s) has a smooth A;^ — > limit). The self-energy 
can be written as 

trP L E(k;l,s)= [ d 4 ph(k,p) exp(-sa(p)) . (17) 

J B.Z. 

The exponential containing the s-dependence arises from the internal ferm- 
ion line with four- momentum p^. All other terms are represented by h(p, k). 

In the limit s ^> 1 one can employ a saddle-point approximation be- 
cause tr Pl S(/c; 1, s) is dominated by the maximum of exp(— a(p)) over the 
Brillouin zone. The exponents a(p) are determined by the tree-level DWF 
action [2, 3]. For ordinary DWF max{exp(— a)} = 0.5 at M = 1. This 
maximum is obtained at the four points P^ = {(ir, 0, 0, 0), (0, ir, 0, 0), . . .}. 
(The values for any < M < 2 were computed in ref. [27].) The s -2 factor 
in eq. (16) arises because the (gaussian) integration is four-dimensional. 

Within tadpole-improved perturbation theory and assuming M has 
been tuned to 1 + SM, the s-dependence of the quark's wave function 
is governed by 5xi{s). Ignoring the pre-exponential factor in eq. (16), we 
thus find to one-loop order the universal result q pt = max{exp(— a)} = 0.5. 
Physically, what the result means is that propagation in the s-direction is 
dominated by the four states belonging to P^. In perturbation theory, the 
wave functions of all quark states have roughly the same s-dependence, 
namely q pt , because they all communicate with the states in P w . 

3.2. NEW FOUR-DIMENSIONAL TERMS 

If we want to improve the fall-off of the quark's wave function, we must look 
for DWF actions where max{exp(— a)} is smaller. We stress that, since it is 
determined by the tree- level DWF action, a(p) is unchanged if one employs 
"fat links" or an improved gauge action. We will first consider here DWF 
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actions where ip x couples to both ip x ±jl and ip x ±2ji, we allow for next- 
nearest neighbors, but only in the same direction. For the s-coordinate we 
retain the same nearest-neighbor coupling as with ordinary DWF. 

The four-dimensional part of the usual DWF action contains the familiar 
Wilson-Dirac operator. Here we replace it by a new operator, whose tree- 
level momentum-space form is 

D 23 = J2^MP») + rW 2- M > ( 18 ) 

/3(p„) = sin(p M ) [1 + c 3 (l - cos(p M ))] , (19) 

W„ = ^(l-cos(^)r, (20) 
ft 

(the standard Wilson-Dirac case corresponds to C3 = 0, n = 1). For M = 1, 
following ref. [2, 3] the exponents are determined by 

2 C ° sh ^ = < 21 ) 

where by convention a > 0. In ref. [13] we analyze the mathematical struc- 
ture of the minima of eq. (21) and explain why taking n = 2 and C3 > 
raises the global minimum of cosh(a) which, in turn, lowers the global max- 
imum of exp(— a). A few examples are given in Table 1. The values in the 
last column should be compared with max{exp(— a)} = 0.5 of ordinary 
DWF. 



TABLE 1. max{exp(— a)} for the operator D23 and for 
various values of C3 and r (see text for the definitions). 
The second column gives the first two terms in the expan- 
sion of fsijp). For each C3 we looked for the r- value where 
max{exp( — a)} is smallest. 



C3 


Hp) 


r 


2min{cosh(a)} 


max{exp(— a)} 




1 

3 
2 
3 
4 
3 
7 
3 


p-sP 3 
p 

P+P 3 


1.46 
1.14 
1.19 
1.45 
1.98 


2.83 
3.40 
4.09 
5.62 
8.06 


0.414 
0.326 
0.261 
0.184 
0.126 



We now digress to discuss how the present work relates to the stan- 
dard "improvement program" (see e.g. the review [28]). In the study of the 
hadron spectrum, only a single parameter (the bare quark mass) in the 
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fermion action needs to be tuned. Once the correct continuum limit has 
been established, attention is focused on eliminating those lattice artifacts 
that vanish most slowly, that is, linearly with the lattice spacing. How- 
ever, in the calculation of weak matrix elements one has to first establish 
the correct continuum limit. This is very problematic with the standard 
fermion methods because, due to the loss of full chiral symmetry, many 
subtraction coefficients must be tuned. Thus, having good chiral properties 
is of higher priority than the removal of any other lattice error. Also, in 
the massless-quark limit, 0(a) lattice artifacts are automatically excluded 
if chiral symmetry is maintained [5]. In that sense, approaching the chiral 
limit using DWF encompasses the standard improvement program as well. 

Coming back to the new DWF action, since the new Wilson term W2 
starts off at order p 4 , the first lattice deviation from a relativistic (tree- level) 
dispersion relation comes only from the kinetic term. This is shown in the 
second column of Table 1. We observe that while increasing C3 from zero 
to 1/3 improves the dispersion relation, the opposite is true for C3 > 1/3. 
Although the error is formally of order a 2 , it might become significant if 
C3 is too large. To gain some idea on the magnitude of the error consider, 
say, p 2 ~ (400 MeV) 2 , which is relevant for kaon physics, on a lattice with 
a -1 ~ 2 GeV. This means a 2 p 2 ~ 1/25. For the last two rows of Table 1, 
the effect is 2% and 4% respectively. 

If next-next-nearest neighbors are also allowed (still only in the same 
direction) one can obtain a smaller max{exp(— a)} while maintaining a 
vanishing j> 3 term. Let 

D 35 = Y, , yi*h(p»)+rW 3 -M, (22) 

where W n is defined in eq. (20) and 
hiPn) = sin(p^) 



1 + g(l - cos(p^)) + c 5 (l - cos(p M )) 2 



(23) 



Some values of max{exp(— a)} are given in Table 2. Considering the last 
row (05 = 50), the deviation from Lorentz covariance is at the level of 
(50/4)(a 2 p 2 ) 2 ~ 2% for the same p 2 . For C5 = 5, the deviation is below 2% 
up to (700 MeV) 2 . 

The values of q pt = max{exp(— a)} in Table 2 are remarkably small. 
Truly, the precise relation between q p i and the non-perturbative q e g is not 
known. Since for ordinary DWF q pt = 0.5 while q e s ~ 0.8, one may hope 
that q pt provides at least an order-of- magnitude estimate of q e g. Even under 
this mild assumption, it seems that the use of D23 or .D35 may allow one to 
reduce N s a lot while maintaining (or even improving) the quality of chiral 
symmetry compared to the ordinary DWF action. 
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TABLE 2. max{exp( — a)} for the operator D 35 and various 
values of C5 and r. The second column gives the first two terms 
in the expansion of /s(p). 



C5 


Mp) 


r 


2min{cosh(a)} 


max{exp(— a)} 


5 

10 
30 
50 


p + ^p 5 
p + 2 wp 5 

P + 7 THP 5 
P + 12^P 5 


2.21 
3.18 
5.70 
7.48 


17.74 
25.44 
45.51 
59.84 


0.0566 
0.0394 
0.0220 
0.0167 



Replacing the four-dimensional part of the DWF action by D23 (-D35) 
approximately doubles (triples) the number of entries in the DWF matrix. 
Therefore one should expect an increase in the cost of simulations by a 
factor of two (three) at fixed N s . Our above considerations suggest that it 
may be possible to compensate for this by working at smaller N s (at least 
in those cases where ordinary DWF require very large N s ). 

In the continuum limit, all the DWF actions considered here support a 
single quark (one chiral zero mode on each boundary) for |1 — M\ < 1. If 
the (generalized) Wilson term W n is employed, there is a four-quark zone 
(corresponding to the doubler states in P n ) at |1 — M + r2 n \ < 1. An 
additional benefit of the new actions is that the four-quark and the single- 
quark zones are separated by a large gap (as a function of M) which we 
expect to persist also in realistic simulations. 

4. CONCLUSION 

While the cost of simulating DWF is much bigger than that of ordinary 
Wilson fermions, existing results show that this is more than compensated 
by the much better quality of chiral symmetry Similarly, employing one 
or both of the modifications proposed here might turn out to have enough 
benefits to outweigh the extra cost of simulating the corresponding DWF 
actions. In particular, it is hoped that these new actions may make it feasi- 
ble to simulate DWF on coarse lattices (a -1 ~ 1 GeV) without loosing the 
good chiral properties. 
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